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i»J ■ deformation quantization, the multiresolution representations and 

the variational approach are the key points. We construct the so- 
lutions of Wigner-like equations via the multiscale expansions in 
the generalized coherent states or high-localized nonlinear eigen- 
j^ ■ modes in the base of the compactly supported wavelets and the 

wavelet packets. We demonstrate the appearance of (stable) lo- 
calized patterns (waveletons) and consider entanglement and de- 
coherence as possible applications. 
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We present the application of the variational- wavelet analysis to the analysis of quantum ensembles in Wigner 
framework. (Naive) deformation quantization, the multiresolution representations and the variational approach 
are the key points. We construct the solutions of Wigner-like equations via the multiscale expansions in the 
generalized coherent states or high-localized nonlinear eigenmodes in the base of the compactly supported wavelets 
and the wavelet packets. We demonstrate the appearance of (stable) localized patterns (waveletons) and consider 
entanglement and decoherence as possible applications. 



1. WIGNER-LIKE EQUATIONS 

In this paper we consider the applications of 
a numerical-analytical technique based on local 
nonlinear harmonic analysis (wavelet analysis, 
generalized coherent states analysis) to the de- 
scription of quantum ensembles. The correspond- 
ing class of individual Hamiltonians has the form 



H(p,q) 



P 
2m 



U( P ,q), 



(1) 



where U{p, q) is an arbitrary polynomial function 
on p, q, and plays the key role in many areas of 
physics [1]. The particular cases, related to some 
physics models, are considered in [2]. Our goals 
are some attempt of classification and the explicit 
numerical-analytical constructions of the existing 
quantum states in the wide class of models. There 
is a hope on the understanding of relation be- 
tween the structure of initial Hamiltonians and 
the possible types of quantum states and the qual- 
itative type of their behaviour. Most important 
in many areas are: localized states, chaotic-like 
or/and entangled patterns, localized (stable) pat- 
terns (waveletons) . Our starting point is the gen- 
eral point of view of a deformation quantization 
approach at least on the Moyal/Weyl/Wigner 
level [1]. In the naive calculations we may use the 



(2) 



simple formal representation for star product: 

In this paper we consider the calculations of the 
Wigner functions W(p,q,t) (WF) corresponding 
to the classical polynomial Hamiltonian H {p, q, t) 
as the solution of the Wigner equation [1] : 

ifi— W = H*W -W*H (3) 

dt v ; 

and related Wigner-like equations. According to 
the Weyl transform, a quantum state (wave func- 
tion or density operator p) corresponds to the 
Wigner function, which is the analogue in some 
sense of classical phase-space distribution [1] . We 
consider the following form of differential equa- 
tions for time-dependent WF, W = W(p, q, t): 



W t = — sin 



-(d H d w 



Hp,W\ 



d"d. 



HW 



(4) 



In quantum statistics the ensemble properties are 
described by the density operator 



P(t) 



Wi|*i(t)><*i(t)|, 






1 (5) 



After Weyl transform we have the following de- 
composition via partial Wigner functions 
Wi(j>, q, t) for the full ensemble Wigner function: 

(6) 



W(p,?,«)=53«xWi(p,g,t) 



where the partial Wigner functions 
W n {q,p,t) = — / d£cxp( - ^p£ 
*;(?-^,*)*n(9+^,*) 



are solutions of proper Wigner equations: 

dw n = P dw n | 

(-i)^(V2) 2£ ^ +1 t/„( g ) a 2£+1 vK„ 

(2*+l)! dq 2 ^ 1 dp u + 1 



(7) 



(8) 



1=0 



The next case describes the important deco- 
hcrcnce process, where we have collective and 
environment subsystems with their own Hilbcrt 
spaces Ti = Ti c <8> T~t e ■ Analysis is based on Weyl 
transform of Lindblad master equation [1]: 



W = {H, W} PB + J2 



n>l 



h 2n (-i) n 

2 2 ™(2n+l)! 



(9) 



d% n+1 U(q)d 2n+1 W(q,p) + 2 1 d pP W + Dd 2 W 

In the next section we consider the variational- 
wavelet approach for the solution of all these 
Wigner-like equations (3), (4), (8), (9) for the case 
of an arbitrary polynomial U(q,p), which corre- 
sponds to a finite number of terms in the series 
expansion in (4), (8), (9) or to proper finite order 
of h. The localized bases/states are the natu- 
ral generalization of standard coherent, squeezed, 
thermal squeezed states [1], which correspond to 
quadratical systems (pure linear dynamics) with 
Gaussian Wigner functions. Representation of 
underlying symmetry group (affinc group in the 
simplest case) on the proper functional space of 
states generate the exact multiscale expansion 
which allows to control contributions to the fi- 
nal result from each scale of resolution from the 
whole underlying infinite scale of spaces. Numeri- 
cal calculations according to methods of part I ex- 
plicitly demonstrate the quantum interference of 
generalized localized states, pattern (entangled- 
like) formation from localized cigenmodes and the 
appearance of (stable) localized patterns (wavele- 
tons). 



2. VARIATIONAL MULTISCALE REP- 
RESENTATIONS 

We obtain our multiscale/multiresolution rep- 
resentations for solutions of Wigner-like equa- 
tions via a variational-wavelet approach. We rep- 
resent the solutions as decomposition into local- 
ized eigenmodes (regarding action of affine group, 
i.e. hidden symmetry of underlying functional 
space of states) related to the hidden underlying 
set of scales [3] : 



W n (t,q,p) = ^Wl(t,q,p) 



(10) 



where value i c corresponds to the coarsest level 
of resolution c or to the internal scale with the 
number c in the full multircsolution decomposi- 
tion of underlying functional space (L 2 , e.g.) cor- 
responding to problem under consideration: V c C 
V c+ i C V c+2 C ... and p = (p 1 ,p 2 ,...), q = 
(q 1 ,q 2 ,...), x t = (pi,qi,---,Pi,qi) are coordi- 
nates in phase space. In the following we may 
consider as fixed as variable numbers of particles. 
The second case corresponds to quantum statisti- 
cal ensemble (via "wignerization" procedure) and 
will be considered in details elsewhere. We intro- 
duce the Fock-like space structure 



H = ($®Hl 



(11) 



for the set of n-partial Wigner functions (states): 

W l = {W*, Wi( Xl ;t), ..., W^xi, ...,x N ;t),.. .} 

where W p (x u ...,x p ;t) e H' } , H° = C, HP = 
L 2 (R ep ) (or any different proper functional spa- 
ce), with the natural Fock space like norm: 



(W, W) = wl + 

J2 [w?(x 1 ,...,x i ;t)'[[ 



(12) 



fie. 



First of all we consider W — W(t) as a func- 
tion of time only, W G L 2 (R), via multires- 
olution decomposition which naturally and effi- 
ciently introduces the infinite sequence of the un- 
derlying hidden scales [3] . We have the contribu- 
tion to the final result from each scale of resolu- 
tion from the whole infinite scale of spaces (11). 



The closed subspace Vj(j e Z) corresponds to 
the level j of resolution, or to the scale j and 
satisfies the following properties: let Dj be the 
orthonormal complement of Vj with respect to 
Vj + i: Vj+i = Vj Dj. Then we have the follow- 
ing decomposition: 



{W(t)}= D j = V c ^D j 



(13) 



— oo<j<oo 



j=o 



in case when V c is the coarsest scale of resolution. 
The subgroup of translations generates a basis 
for the fixed scale number: spa,n k£Z {2 J ' 2 ^(2 J t — 
k)} = Dj. The whole basis is generated by action 
of the full affine group: 



: {2 3 ' 2 ^{2H-k)} 



spa,n keZ j eZ ] 
8pan fciiez {*j-, fc } = {W{t)} 



(14) 



One of the key points (the so called Fast Wavelet 
Transform, FWT) of wavelet analysis approach 
demonstrates that for a large class of operators 
wavelets are good approximation for true eigen- 
vectors and the corresponding matrices are al- 
most diagonal. FWT gives the maximum sparse 
form of operators under consideration [3]. So, 
let us denote our (integral/differential) opera- 
tor from equations under consideration as T 
(L 2 (R n ) -» L 2 {R n )) and its kernel as K. We 
have the following representation: 



<Tf,g>-- 



K{x,y)f(y)g(x)dxdy (15) 



In case when / and g are wavelets ipj.k = 
2 3 ' 2 (p(2 3 x — k) (15) provides the standard rep- 
resentation for operator T. Let us consider mul- 
tiresolution representation . . . C Vi C V\ C Vq C 
V-i C V-2 ■ ■ ■■ The basis in each Vj is ipj t k(x), 
where indices k, j represent translations and scal- 
ing respectively. Let Pj : L 2 (R n ) — ► Vj (j E Z) 
be projection operators on the subspace Vj cor- 
responding to level j of resolution: (Pjf)(x) = 



E/t < f>Vj,k > v 3 A x )- Let Qj 



Pj-i - Pj 



be the projection operator on the subspace Dj 
(Vj-i = Vj ® Dj) then we have the following 
representation of operator T which takes into ac- 
count contributions from each level of resolution 
from different scales starting with the coarsest 



and ending to the finest scales [3]: 

t = Y^iQ^Qi + Qo TP J + p 3 T Qi)- 

3&Z 



(16) 



We remember that this is a result of presence of 
affine group inside this construction. The non- 
standard form of operator representation [3] is a 
representation of operator T as a chain of triples 
T = {Aj,Bj 7 Tj}j e z, acting on the subspaces Vj 



and Dj: Aj 



Di 



Dj,Bj 



Vj 



Dj,Tj,Dj 



Vj, where operators {Aj, Bj,Tj}jt=z arc defined 

as J\j ^ LJ j ± (^/ j , £}j ^ LJ j ± ±j, 1 j -— r j ± C^t. 

The operator T admits a recursive definition via 



Ti 



A 






B 



i+i 



T 



j+i 



where T,- 



PjTPj and Tj acts on Vj : Vj 



(17) 



Vj- 



So, it is possible to provide the following "sparse" 
action of operator Tj on sufficiently smooth func- 
tion /: 



kez \ e 



'j,k-e ¥>j,k(x), (18) 



in the wavelet basis ipj^ix) = 2 ^ 2 cp(2 3 x — k) 
where 

fj,k-i = 2^' /2 / f(x)ip(2- j x -k + £)dx (19) 



are wavelet coefficients and rg are roots of some 
additional linear system of equations related to 
the "type of localization" [3] . So, we have simple 
linear parametrization of matrix representation of 
our operators in localized wavelet bases and of the 
action of this operator on arbitrary vector/state 
in proper functional space. The variational ap- 
proach [2] reduces the initial problem to the prob- 
lem of solution of functional equations at the first 
stage and some algebraical problems at the sec- 
ond one. So, the solution is parametrized by the 
solutions of two sets of reduced algebraical prob- 
lems, one is linear or nonlinear (depending on the 
structure of the generic operator L [2]) and the 
rest are linear problems related to the computa- 
tion of the coefficients of the Galerkin-like alge- 
braic equations [3]. As a result the solution of 
the equations from Sec. 1 has the following mul- 
tiscale decomposition via nonlinear high-localized 



eigenmodes 

W(t,xi,x 2 ,...)= Yl aijU i ®V j (t,xi,...), 

y j{t) = v J,sio W(t) + J2 vifat), ^ ~ 2 l , (20) 



1>N 



u i (x s ) = u i f° w (x s y 



/ , Um\™m X s)> k m ~ 2 



i>A/ 



which corresponds to the full multiresolution ex- 
pansion in all underlying time/space scales. The 
formulae (20) give the expansion into a slow part 
and fast oscillating parts for arbitrary N, M. We 
may move from the coarse scales of resolution to 
the finest ones for obtaining more detailed infor- 
mation about the dynamical process. In this way 
one obtains contributions to the full solution from 
each scale of resolution or each time/space scale 
or from each nonlinear eigenmode. It should be 
noted that such representations give the best pos- 
sible localization properties in the corresponding 
(phase)space/time coordinates. Formulae (20) do 
not use perturbation techniques or linearization 
procedures. Numerical calculations are based on 
compactly supported wavelets and wavelet pack- 
ets and on evaluation of the accuracy on the level 
N of the corresponding cut-off of the full sys- 
tem regarding norm (12): \\W N+1 - W N \\ < e. 
Our (nonlinear) eigenmodes are more realistic 




Figure 1. Localized pattern- like (waveleton) 
Wigner function. 




Figure 2. Entangled- like Wigner function. 



for the modelling of nonlinear classical/quantum 
dynamical process than the corresponding linear 
gaussian-like coherent states. Here we mention 
only the best convergence properties of the ex- 
pansions based on wavelet packets, which real- 
ize the minimal Shannon entropy property and 
the exponential control of convergence of expan- 
sions like (20) based on the norm (12). Figures 
1, 2 present the solutions, constructed from the 
first 6 eigenmodes (6 levels in formula (20)), and 
demonstrate the stable localized pattern forma- 
tion (waveleton) and complex entangled-like be- 
haviour. Fig. 1 corresponds to (possible) result of 
superselection (einselection) [1] after decoherence 
process started from entangled pattern demon- 
strated on Fig. 2. It should be noted that we can 
control the type of behaviour on the level of the 
reduced algebraical variational system [2] . 
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